Electronic Structure and Bulk Spin Valve Behavior in Ca 3 Ru 2 07 
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We report density functional calculations of the magnetic properties and Fermiology of Ca3Ru207. 
The ground state consists of ferromagnetic bilayers, stacked antiferromagnetically. The bilayers are 
almost but not exactly half-metallic. In the ferromagnetic state opposite spin polarizations are 
found for in-plane and out-of-plane transport. Relatively high out of plane conductivity is found 
for the majority spin, which is relatively weakly conductive in-plane. In the ground state in-plane 
quantities are essentially the same, but the out of plane transport is strongly reduced. 



The perovskite based ruthenates, £> 1+ „Ru„03 n+1 , 
Z?=Sr,Ca show a remarkable range of electronic and mag- 
netic properties, even though they are all based on Ru 4+ 
in octahedral environments with corner sharing topolo- 
gies. This includes robust itinerant ferromagnetism 
(SrRu03)i£, paramagnetic "bad" metals (CaRuOs) 3 , 
unconventional superconductivity (S^RuC^) 4 , an anti- 
ferromagnetic (AFM) Mott insulator (Ca2Ru04)£, and 
metamagnetic quantum critical behavior (Sr 3 Ru207) 6 . 
This reflects exceptionally strong dependence of elec- 
tronic and magnetic properties on lattice degrees of free- 
dom, also seen in band structure studies 7 * 8 **!* 1 ^. 

The bilayer n — 2 compounds are of particular interest 
due to their borderline properties. Sr3Rii207 bridges the 
metallic ferromagnetic (FM) n — oo compound SrRuC>3 
and the paramagnetic Fermi liquid Sr2Ru04 and has a 
metamagnetic quantum critical poinfeS associated with 
borderline metallic ferromagnetism. Ca3Ru2C>7 is inter- 
mediate between the bad metal CaRuC>3 and the Mott 
insulator Ca2RuC>4, and so may be a useful window 
into the physics of clean materials near a metal insu- 
lator transition. Signatures of the borderline physics in- 
clude strong sensitivity of measured properties to minor 
sample variations, observation of metallic like proper- 
ties such as quantum oscillations and finite linear specific 
heat coefficient in non-conducting material, a large lat- 
tice anomaly at the magnetic ordering temperature, and 
low field metamagnetic transitions with strong transport 
signaturea 11 . 12 ! 13 . 14 ! 15 . 16 ! 17 . 18 ! 19 ! 20 ! 21 . 

Here, we report detailed electronic structure calcula- 
tions for Ca3Ru2C>7. These give insight into a number 
of key observations - magnetic structure, lattice cou- 
pling, and transport. They were done in the local spin 
density approximation (LSDA) using the general poten- 
tial linearized augmented planewave (LAPW) method 
with local orbitalsSSi 2 ^. We used the experimental low 
temperature (8K) non-centrosymmetric Bb1\m crystal 
structure^. This is a centered structure, with two for- 
mula units, and a c(2x2) in-plane cell doubling due to 
octahedral tilts and rotation. The strong octahedral tilt- 
ing differentiates this compound from Sr3Ru2 07. Sig- 
nificantly, if insulating, this non-centrosymmetric space- 
group would have a ferroelectric polarization along b. 
This would be of interest in the context of magneto- 
electrics if a related well insulating, but magnetic, ma- 



terial can be made. The calculations were highly con- 
verged with respect to basis set and zone samplings 2 ^*^. 
Approximately 2650 LAPW basis functions were used, 
and twice that number for the doubled cell. The Fermi 
surfaces are based on an interpolation from 494 first prin- 
ciples points in the irreducible wedge of the zone, except 
for the doubled AFM cell, where 194 points were used. 

Ca3Ru2C>7 shows an AFM ordering at T/v=56K and a 
transition from a metal to a low temperature poorly con- 
ductive or insulating phase at 48K. 11 Metamagnetic tran- 
sitions to an effectively ferromagnetic phase, with large 
changes in conductance occur at relatively low field, de- 
pending on the field direction while magnetization mea- 
surements imply that the RuC>2 layers themselves are fer- 
romagnetically ordered, with an AFM stacking 1 ^. Three 
possible AFM stackings have been suggested - ferro- 
magnetic bi-layers stacked antiferromagnetically, or bi- 
layers that arc internally antiferromagnetically aligned, 
stacked either in a FM or an AFM fashion, i.e. (1) 
-UU - DD - UU-, (2) -UD - UD - UD- or (3) 
—UD — DU — UD— along c. Recent neutron scattering 
experiments favor ordering (1) over ordering (2), based 
on observation of half order diffraction peakst^i Consid- 
ering the layered crystal structure and highly 2D elec- 
tronic structure (see below) , magnetic interactions within 
a bilayer must be very much stronger than inter-bilayer 
couplings. Calculated energetics are given in Table [I] 
The non- spin-polarized state (P) is unstable against fer- 
romagnetism by 70 meV/f.u., but state (2), is favored 
over P by only 19 meV/f.u. In plane antiferromagnetism 
is also disfavored relative to ferromagnetism. Thus AFM 
in the bilayers is strongly disfavored, and the ground 
state must be (1) and not (3). In fact, (1), consisting 
of ferromagnetic bi-layers stacked antiferromagnetically 
is an additional 5 meV/f.u. lower than the FM state, so 
the inter-bilayer coupling is weakly AFM. We did struc- 
tural relaxations of all atomic positions in the Bb2\m 
cell. This yields an energy gain of only 16 meV (per 
12 atom formula unit) for FM ordering, supporting the 
experimental crystal structure, and a larger 27 meV for 
the non-spin-polarized case, the difference implying sig- 
nificant magnetoelastic coupling. For magnetism with 
itinerant character, this can provide an explanation of 
the structural anomalies at the ordering temperature. 

Since the inter-bilayer coupling is weak, it is useful to 
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TABLE I: Calculated magnetic energies. AF1 and AF2 are 
FM layers stacked -UU-DD-UU- and -UD-UD-UD-, respec- 
tively. AFP is in-plane c(2x2) AFM. 



ordering 


P 


FM 


AFP 


AF1 


AF2 


E(meV/f.u.) 





-72 


-52 


-77 


-19 



start with the FM state. The calculations show that the 
mechanism of moment formation and the major sheets 
of Fermi surface are essentially the same as in the AFM 
ground state, as are the density of states and in-plane 
Fermi velocities. However, for the AFM case, the major- 
ity and minority spin lie on top of each other due to the 
alternating spin up and spin down bi-layers. 

The spin magnetization with FM ordering is 1.934 
/is/Ru, of which only 1.23 /is is within the Ru LAPW 
sphere (radius 2.05 oq). The remaining ~ 1/3 is O in 
character, comparable to other FM ruthenates 8,9 . The 
origin is the strong Ru-0 hybridization in these 4d com- 
pounds, combined with Hund's coupling on O. The elec- 
tronic density of states (DOS) and projections onto the 
Ru and O LAPW spheres are shown in Figs. ^ an d ID 
The optical spectrum is in Fig. [3] and the Fermi surface 
in Fig. 0] The low energy part of the optical spectrum 
below ~ 1.5 eV, derives from d-d transitions, while the 
in-plane peaks at ~ 1.9 eV and ~ 2.5 eV are respectively 
of majority and minority spin charge transfer p-d char- 
acter. There is little in-plane anisotropy except at low 
energy (not shown) . Comparison with experiment would 
be useful in testing the LSDA electronic structure. 

Moment formation via an itinerant Stoner instability 
depending on Ru-0 hybridization, similar to SrRu03 is 
clearly seen££. Thus the O polarization plays a key role. 
The importance of the Hund's coupling on O is also evi- 
dent from the exchange splitting of the O bands as seen in 
the DOS. The spin moments within the O LAPW spheres 
(radius 1.65 ao) are 0.08 hb (apical O), 0.13 /ib and 0.14 
\xb (two inequivalent in-plane O) and 0.21 fis (interlayer 
bridging O). These are reflected in the various O contri- 
butions to the DOS near Ep (Fig. |2J. The large bridging 
O contribution, which also reflects the bonding of the bi- 
layer units, explains the strong interlayer FM coupling. 

Our results and the experimental observation of rela- 
tively low field metamagnetism show that the basic mag- 
netic structure is ferromagnetic bilayers, with a very weak 
inter-bilayer AFM coupling. To understand the elec- 
tronic structure we begin with the FM bilayers. 

Within an ionic model, each Ru has 4 ti g electrons, 
which partially fill the t2 g manifold. Because of band 
narrowing due to the octahedral tilts and rotation in 
Ca3Ru2 07, this manifold is separated from the higher 
lying e g manifold by a gap. If fully polarized, this would 
yield a spin moment of 2 /x^/Ru, close to but larger than 
what we find. The FM DOS shows that the bilayers 
are almost, but not exactly half metallic. This is simi- 
lar to the colossal magnetoresistive (CMR) manganites 
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FIG. 1: (color online) Electronic density of states and pro- 
jection onto the Ru LAPW spheres, radius 2.05 ao) of FM 
CaaR^Cv on a per formula unit basis. Majority spin is shown 
above the axis, and minority below. 
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FIG. 2: (color online) Density of states projection onto the 
O LAPW spheres (radius 1.65 ao) on a per atom basis. 01 
is the apical O in the CaO rocksalt layers, 02 is the bridging 
apical O joining the Ru02 bilayers and 03 and 04 are the 
plane O in the Ru02 layers. 



such as (La,Ca)MnOg2Si21. In those materials, there is 
a finite density of states in both spin channels, but Ef 
falls near the minority spin band edge, and those states 
are Anderson localized. In FM Ca3Ru20y, Ep falls very 
near the majority spin band edge. However, because it 
is a clean material, with no alloy scattering, there is no 
mechanism for Anderson localization at low temperature. 
Thus, in contrast to the CMR materials, as the tempera- 
ture is lowered in FM Ca 3 Ru207 the majority spin chan- 
nel will cross-over from localized to metallic, yielding an 
anisotropic partially spin polarized metal. 

In a layered structure, the ti g manifold derives from 
two sets of bands, d xy which is 2D, and d xz and d yz , 
which are ID. For a bilayer, these split into symmetric 
and antisymmetric combinations due to interlayer hop- 
ping. Since the d xy is directed in-plane, its interlayer 



TABLE II: Fermi surface parameters of Ca3Ru207 for the P, 
FM and AFM ground state AF1. N(E F ) is given in eV _1 per 
spin per f.u. Fermi velocities, < wf > 1//2 are in cm/s. 





N(Er) 






Vc 


p 


5.23 


0.78xl0 7 


0.71xl0 7 


0.23xl0 7 


F(majority) 


0.95 


0.64xl0 7 


0.71xl0 7 


0.58xl0 7 


F(minority) 


4.41 


1.17xlO r 


1.48xl0 7 


0.22xl0 7 


AF1 


2.76 


l.OxlO 7 


1.4xl0 7 


O.lOxlO 7 



coupling is expected to be smaller than for the d xz and 
d yz bands. Conversely, its in-plane band width may be 
expected to be larger, following the number of hopping 
paths. This qualitatively is what we find. 

In Boltzmann theory with the constant scattering 
time approximation, the conductivity varies as a x ~ 
N(Ef)v Fx t, where r is a scattering time. Both the 
paramagnetic and FM electronic structures are quite 
anisotropic near Ep as expected from the layered crystal 
structure. This is also apparent from the Fermi surfaces. 
The exception is the majority spin in the FM case, which 
has a low density of states and is 3D in character. This 
Fermi surface has mixed character due to tilting, but is 
largely from the bilayer antisymmetric combination of 
d xy orbitals, Hopping through the CaO rocksalt layers 
is geometrically allowed by octahedral tilts in this struc- 
ture, so that even though the Fermi surface has mostly 
d xy character it involves enough apical O character to 
produce c-axis dispersion. 

Turning to the minority spin channel, in addition to 
the highly 2D character, there is also significant in plane 
anisotropy. In the constant scattering time approxima- 
tion, b direction conductivity is ^50% higher than along 
a. This means that sample dependent orthorhombic 
twinning could significantly affect the measured trans- 
port and should be considered in interpreting experi- 
ments. In plane conduction is dominated by the minority 
channel and is highly spin polarized 2 ^, P^ = -0.88 and 

(2) 

P h = -0.91. On the other hand, for the FM ordering, 
while c-axis conductivity is a factor of 20 lower, it has 
a small opposite polarization, Pc — +0.19. The actual 

(2) 

value of Pc may actually be somewhat larger, because 
the multi-sheet minority spin Fermi surface would pro- 
vide more non-spin-flip scattering channels than the sim- 
ple majority surface, and so the minority spin scattering 
time may be shorter than the majority spin. In any case, 
to our knowledge, this type of direction dependent sign 
of the polarization is unique to FM Ca3Ru 2 0y. While 
here it is a scientific curiosity, if found in materials that 
have true ferromagnetic ground states and large \P\ for 
both directions, it could be of practical importance for 
sources of spin polarized electrons with spin depending 
on the contact geometry. 

The majority spin Fermi surface consists of a single 
highly corrugated hole cylinder of mostly d xy charac- 
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FIG. 3: Optical spectrum (e 2 ) for FM Ca 3 Ru 2 7 . 




FIG. 4: (color online) Majority (top) and minority (bottom) 
basal plane extended zone FM Fermi surfaces. The actual 
zone is half the area and rotated 45°. The plot shows a color 
map ranging from 1 mRy below E F (blue) to 1 mRy above 
(red). The width is inversely proportional to the velocity and 
electron (hole) surfaces are blue (red) on the interior. The 
center is a Z point; alternating corners are V and Z points 
due to the centered lattice. 
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ter around the zone center, which almost pinches off at 
k z = 0. The minority spin has several sheets of Fermi 
surface. The highest velocity parts, which contribute 
most to the in-plane conductivity, are from the symmet- 
ric and antisymmetric ID d xz and d yz sections. These 
reconnect to form a square cylindrical T centered section 
and an open section, that forms from a square cylinder 
but opened at the zone boundary at (1/2,0), which is 
(1/4,1/4) in the extended zone of Fig. Along the 
other 6-direction zone boundary there is a complex set 
of Fermi surfaces. These derive from the 2D d xy sheets 
which kiss the zone boundary at (0,1/2). These do not 
appear along the a direction because the d xy sheet is 
gapped away by mixing with the ID d xz and d yz bands 
due to the octahedral tilts, which are in the a-direction 
(i.e. around the 6-axis). The repulsion from the d xz /d yz 
Fermi surfaces is strong enough to fold this sheet back 
in the region around (1/2,1/2) leading to the complex 
structure seen around this point. The presence of the 
d xy derived sections towards the 6-direction zone bound- 
ary, but not the a-direction, is the reason for the large 
in-plane anisotropy of a. Besides the major sheets there 
are two smaller concentric nearly circular cylindrical elec- 
tron sections around T, which are also of substantial d xy 
character. Significantly, these two minority spin sheets 
have the same center and are intermediate in size be- 
tween the maximum and minimum size of the majority 
spin Fermi surface as a function of k z . However, they are 
mis-matched in that the majority surface contains holes 
and the minority contains electrons. 



Turning to the AFM ground state, as mentioned, one 
might expect the electronic structure of the bi-layers to 
be essentially unchanged, due to the weak inter-bilayer 
coupling, and this is confirmed by the calculations. The 
only notable change is a disruption of the c-axis trans- 
port. There are no significant changes^ in N(Ep) or the 
in-plane N(Ep)v F . However, the spin averaged c-axis 
N(Ep)vp is strongly reduced by a factor of 10 relative 
to the FM case. This is strong spin-valve physics. More- 
over, there is another factor. For non-spin- flip scatter- 
ing, which ordinarily dominates, there are more scatter- 
ing channels for the minority spin due to its larger multi- 
sheet Fermi surface. Thus with AFM stacking, there may 
be an important change in r, which would further lower 
the c-axis conductivity, depending on the relative impor- 
tance of spin-flip and non-spin-flip scattering, which in 
turn would depend on details of the sample and temper- 
ature. 

In summery, we report calculations of the electronic 
and magnetic properties of Ca3Ru207, which elucidate 
the ground state ordering and provide a framework 
for the strong connections between magnetic order and 
transport properties, particularly at the metamagnetic 
transition. 
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Department of Energy. 
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